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Abstract. 

Supersymmetric models provide very interesting scenarios to account for the dark matter of the Universe. In this talk 
we discuss scenarios with gravitino dark matter in /^-parity breaking vacua, which not only reproduce very naturally the 
observed dark matter relic density, but also lead to a thermal history of the Universe consistent with the observed abundances 
of primordial elements and the observed matter-antimatter asymmetry. In this class of scenarios the dark matter gravitinos are 
no longer stable, but decay with very long lifetimes into Standard Model particles, thus opening the possibility of their indirect 
detection. We have computed the expected contribution from gravitino decay to the primary cosmic rays and we have found 
that a gravitino with a mass of ra3/2 ~ 150 GeV and a lifetime of ~ 10^* s could simultaneously explain the EGRET 
anomaly in the extragalactic gamma-ray background and the HEAT excess in the positron fraction. 
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INTRODUCTION 

A series of observations have provided in recent years 
compelling evidences for the existence of dark matter in 
the Universe These observations have revealed that 
the dark matter particle has to be weakly interacting with 
the ordinary matter, long lived and slow moving ("cold"). 
Among the most interesting candidates proposed stands 
the gravitino |2] which is abundantly produced by ther- 
mal scatterings in the very early Universe. If the gravitino 
is the lightest supersymmetric particle (LSP), a fraction 
of their initial population will survive until today with 
a relic density which is calculable in terms of very few 
parameters, the result being ^ 

^'^ VlO'^GeVyV '"3/2 / VlTeV/ ' 

(1) 

where Tr is the reheating temperature of the Universe, 
m3/2 is the gravitino mass and is the gluino mass. 
In predicting the relic abundance of gravitinos, the main 
uncertainty arises from our ignorance of the thermal 
history of the Universe before Big Bang nucleosynthesis 
(BBN) and in particular of the reheating temperature 
after inflation. However, we have strong indications that 
the Universe was very hot after inflation. Namely, the 
discovery of neutrino masses provided strong support to 
leptogenesis as the explanation for the observed baryon 
asymmetry of the Universe |i4J. This mechanism can 
reproduce the observed baryon asymmetry very naturally 
if the reheating temperature of the Universe was above 
10^ GeV [5]. Therefore, the abundance of gravitinos 
can reproduce the dark matter relic density inferred by 



WMAP for the ACDM model, ilcDM/j^ - 0.1 R], for 
natural values of the input parameters in Eq. ([T]). 

Being capable of reproducing the correct relic density 
is a necessary requirement for any dark matter candi- 
date, but not the only one. Namely, the physical model 
accounting for the dark matter should not spoil the suc- 
cessful predictions of the standard cosmology. However, 
in scenarios with gravitino dark matter, when 7? -parity is 
conserved, the next-to-LSP (NLSP) is typically present 
during or after Big Bang nucleosynthesis, jeopardizing 
the successful predictions of the standard nucleosynthe- 
sis scenario. This is in fact the case for the most likely 
candidates for the NLSP: the lightest neutralino and the 
right-handed stau. More precisely, when the NLSP is the 
neutralino, the hadrons produced in the neutralino decays 
typically dissociate the primordial elements [7], yielding 
abundances in conflict with observations. On the other 
hand, when the NLSP is a charged particle, , the for- 
mation of the bound state i^W&X^^ catalyzes the produc- 
tion of ''Li 1 8] leading to an overproduction of ^Li by a 
factor 300- 600 [9J. 

Different solutions have been proposed to this prob- 
lem. For instance, in some specific supersymmetric mod- 
els the NLSP can be a sneutrino |1^ or a stop fllil , 
whose late decays do not substantially affect the predic- 
tions of Big Bang nucleosynthesis. Other solutions are to 
assume a large left-right mixing of the stau mass eigen- 
states lfl2ll or to assume some amount of entropy pro- 
duction after NLSP decoupling, which dilutes the NLSP 
abundance ifTsll . Our proposed solution consists in intro- 
ducing a small amount of 7? -parity violation, so that the 
NLSP decays into two Standard Model particles before 



the onset of Big Bang nucleosynthesis, thus avoiding the 
BBN constraints altogether lfl411 . 



GRAVITINO DECAY 

When /^-parity is not exactly conserved, the gravitino 
LSP is no longer stable. Nevertheless, the gravitino decay 
rate is doubly suppressed by the Planck mass and by the 
small /^-parity violation 1 15]. Therefore, for the range of 
/^-parity violating couplings favored by cosmology, the 
gravitino lifetime ranges between 10^"' and 10^^ s for 
"^3/2 = 100 GeV, which exceeds the age of the Universe 
by many orders of magnitude. Hence, even though the 
gravitino is not absolutely stable, it is stable enough to 
constitute a viable candidate for the dark matter of the 
Universe, while preserving the attractive features of the 
standard Big Bang nucleosynthesis scenario and thermal 
leptogenesis. 

Interestingly, gravitinos could be decaying at a rate 
sufficiently large to allow the detection of the decay 
products in experiments, thus opening the possibility of 
the indirect detection of the elusive gravitino dark matter. 

In this talk we will discuss the possibility that grav- 
itinos are heavier than the gauge bosons. If this is the 
case, gravitinos decay mainly into three different chan- 
nels iSlnlH: Wi/i TV, W^t^, tPv. The fragmen- 
tation of the and the Z" eventually produces a flux 
of stable particles, which we have simulated with the 
event generator PYTHIA 6.4 JTl] . On the other hand, the 
branching ratios of the different decay channels are also 
calculable in the framework of supergravity with broken 
/^-parity, yielding a result that depends mainly on the 
gravitino mass (under the popular assumption of gaug- 
ino mass unification at the Grand Unification scale). 

Dark matter gravitinos populate the halo with a dis- 
tribution that follows a density profile p(r), where r 
denotes the position with respect to the center of the 
Galaxy. We will adopt for our numerical analysis the 
Navarro-Frenk- White density profile ll20ll : 



Po 



(rA,)[l + (r/r,)]2 



(2) 



with Po ~ 0.26GeV/cm^ and ~ 20 kpc, although our 
conclusions are not very sensitive to the choice of the 
density profile. Gravitinos at r eventually decay with 
lifetime T3/2 producing photons, positrons, antiprotons 
and neutrinos at a rate per unit energy and unit volume 
given by 



(3) 



where x = 7, e+ , , v and dNx/ dE is the energy spectrum 
of the particle x produced in the decay. Remarkably, 



the source function depends essentially on two unknown 
parameters, namely the gravitino mass and the gravitino 
lifetime, yielding a fairly predictive scenario. 

We will present in this talk the results for the 
gamma-ray flux II6I1 and for the positron and antipro- 
ton fluxes iITtIi . The results for the neutrino flux will be 
presented elsewhere ll2lll . The flux at Earth of the dif- 
ferent particle species are constrained by a series of ex- 
periments. EGRET measured gamma-rays in the energy 
range between 30 MeV to 100 GeV. After subtracting 
the galactic foreground emission, the residual flux was 
found to be roughly isotropic and thus attributed to ex- 
tragalactic sources. An improved analysis of the galac- 
tic foreground by Strong et al. fl^, optimized in order 
to reproduce the galactic emission, shows a power law 
behavior between 50 Me V-2 GeV, but a clear excess be- 
tween 2-10 GeV, roughly the same energy range where 
one would expect a signal from gravitino decay. On the 
other hand, the flux of positrons has been measured by 
a series of experiments, in particular by HEAT |23]. Al- 
though the measurements still suffer from large uncer- 
tainties, it is intriguing that all the experiments seem to 
point to an excess of positrons at energies larger than 7 
GeV, which is again the energy range where a contribu- 
tion to the flux from gravitino decay is expected. Lastly, 
the measurements of the antiproton flux by BESS ll24ll 
and other experiments do not show any deviation from 
the predictions by conventional astrophysical models of 
spallation of cosmic rays on the Milky Way disk. Future 
experiments such as the FGST (formerly GLAST) ll25ll . 
measuring the gamma-ray flux, and PAMELA 1^2^, mea- 
suring antimatter fluxes, will provide in the near future 
very precise measurements of the spectra of cosmic rays 
which will constitute decisive tests of the decaying dark 
matter scenario. 



GAMMA-RAY FLUX 

The gamma-ray flux from gravitino decay has two com- 
ponents. The decay of gravitinos at cosmological dis- 
tances will be detected at Earth as a perfectly isotropic 
diffuse gamma-ray background with a red-shifted energy 
spectrum. A second source of gamma-rays is the decay of 
gravitinos in the Milky Way halo. For typical halo mod- 
els, we find that the halo component dominates over the 
cosmological one lll4ll27ll . yielding a slightly anisotropic 
gamma-ray flux. 

The different contributions to the total gamma-ray 
flux from gravitino decay are shown in Fig. [T] for a 
mass of m3/2 ^ 150 GeV and a lifetime of T3/2 1-3 x 
10^6 s d]. To compare our results with the EGRET 
data ll22ll . also shown in the figure, we have averaged 
the halo signal over the whole sky excluding a band of 




FIGURE 1. Contributions to the total gamma-ray flux for 
'"3/2 = 150 GeV and T3/2 — 1-3 x 10^^ s compared to the 
EGRET data. In dotted lines we show the photon flux from 
the fragmentation of the Z boson, in dashed lines from the 
fragmentation of the W boson, and in dot-dashed lines from 
the two body decay 11/3/2 jv. The background is shown as 
a long dashed line, while the total flux received is shown as a 
thick solid line. 



±10° around the Galactic disk, and we have used an en- 
ergy resolution of 15%, as quoted by the EGRET col- 
laboration in this energy range. The energy resolution 
of the detector is particularly important to determine the 
width and the height of the monochromatic line stem- 
ming from the two body decay \//3/2 yv. The three 
contributions are dominated by the halo component, the 
extragalactic component being smaller by a factor of 2- 
3. Finally, we have adopted an energy spectrum for the 
extragalactic background described by the power law 
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(cm str s GeV) , in or- 



der to provide a good agreement of the total flux received 
with the data. 

The predicted energy spectrum shows two qualita- 
tively different features. At energies between 1-10 GeV, 
we expect a continuous spectrum of photons coming 
from the fragmentation of the gauge bosons. As a result, 
the predicted spectrum shows a departure from the power 
law in this energy range that might be part of the appar- 
ent excess infeiTed from the EGRET data by Strong et 
al. ll22ll . In addition to the continuous component, the en- 
ergy spectrum shows a relatively intense monochromatic 
line at higher energies arising from the decay channel 

V/3/2 ^ yv. 

Even though the gamma-ray flux from gravitino decay 
is expected to be anisotropic, we find that with these 
parameters the flux resembles an isotropic extragalactic 
flux in EGRET. Namely, in the energy range 0.1-10 
GeV, the anisotropy between the Inner Galaxy region 
{\b\ > 10°, 270° < I < 90°) and the Outer Galaxy region 
i\b\ > 10°, 90° < / < 270°) is just a 6%, wefl compatible 
with the EGRET data iH. 



ANTIMATTER FLUX 

Antimatter propagation in the Milky Way is commonly 
described by a stationary two-zone diffusion model with 
cylindrical boundary conditions |28]. Neglecting reac- 
celeration effects and non-annihilating interactions of 
antimatter in the Galactic disk, the propagation can be 
described by only a few parameters, which can be de- 
termined from the flux measurements of other cosmic 
ray species, mainly from the Boron to Carbon (B/C) ra- 
tio ll29tl . 

The solution of the transport equation at the Solar 
System can be formally expressed by the convolution 



fiT) 
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dT'G{Tj' 



dN{T') 



(4) 



'W3/2T3/2 Jo ~ ' dT' 

where T is the kinetic energy and T^^ax = 'W3/2 for the 
case of the positrons while Zrnax = 'W3/2 - m,, for the an- 
tiprotons. The solution is thus factorized into two parts. 
The first part, given by the Green's function G{T, T'), en- 
codes all of the information about the astrophysics (such 
as the details of the halo profile and the complicated 
propagation of antiparticles in the Galaxy) and is univer- 
sal for any decaying dark matter candidate. The remain- 
ing part depends exclusively on the nature and properties 
of the decaying dark matter candidate, namely the mass, 
the Ufetime and the energy spectrum of antiparticles pro- 
duced in the decay. Finally, the flux of primary antiparti- 
cles at the Solar System from dark matter decay is given 
by: 



An 



(5) 



where v is the velocity of the antimatter particle. 

Clearly, if gravitino decay is the explanation for the 
extragalactic EGRET anomaly, our predicted positron 
and antiproton fluxes should not exceed the measured 
ones. In the scenario we are considering the gravitino 
mass and lifetime are constrained by requiring a qual- 
itatively good agreement of the predicted extragalactic 
gamma-ray flux with the EGRET data: ;7i3/2 = 150 GeV 
and T3/2 = 1.3 X lO^^s iH]. On the other hand, the en- 
ergy spectrum of antiparticles, dN /dT, is determined by 
the physics of fragmentation. Therefore, the main uncer- 
tainty in the computation of the antimatter fluxes stems 
from the determination of the Green's function, i.e. from 
the uncertainties in the propagation parameters and the 
halo profile. We have found in our analysis that the un- 
certainties in the precise shape of the halo profile are not 
crucial for the determination of the primary antimatter 
fluxes. On the other hand, the uncertainties in the propa- 
gation parameters can substantially change the predic- 
tions for the antimatter fluxes, even by two orders of 
magnitude for the antiproton flux. 

Positrons and antiprotons have different properties and 
their respective transport equations can be approximated 



TABLE 1. Coefficients of the interpo- 
lating function Eq. l|6j for tlie positron 
Green's function, assuming a NFW halo 
profile and for the M2, MED and Ml 
propagation models defined in |30J. 



model 


a 


b 


M2 


-0.9716 


-10.012 


MED 


-1.0203 


-1.4493 


Ml 


-0.9809 


-1.1456 



by different limits of the transport equation. Let us dis- 
cuss each case separately. 

Positron flux 

In the case of positron propagation, galactic conven- 
tion and annihilations in the disk can be neglected in the 
transport equation. We have solved the transport equation 
and computed the Green's function for three propagation 
models, denoted as M2, MED and Ml, which are consis- 
tent with the B/C ratio and which yield, respectively, the 
minimum, median and maximal positron flux |30]. The 
Green's function from dark matter decay can be well ap- 
proximated by the following interpolating function llTll : 
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0(r'-r)cm"^s , 



(6) 

where T and T' are expressed in units of GeV. The co- 
efficients a and b for each propagation model can be 
found in Table [T] for the NFW profile. This approxima- 
tion works better than a 15-20% over the whole range of 
energies. It should be stressed that this parametrization 
of the Green's function is vaUd for any decaying dark 
matter particle, not just for the gravitino. 

With the previous parametrization of the Green's func- 
tion, it is straightforward to compute the interstellar 
positron flux using Eqs. ^ and Eq. (|5]l. The total 
positron flux received at Earth receives, in addition to the 
primary flux from gravitino decay, a secondary compo- 
nent originating in the collision of primary protons and 
other nuclei on the interstellar medium, which constitutes 
the background to our signal. 

To compare our predicted flux with the observations 
we choose to show the positron fraction, defined as the 
total positron flux divided by the total electron plus 
positron fluxes: 



PF(r) 



(7) 

where following IBU 13211 we have left the normalization 
of the primary electron flux, k, as a free parameter to be 
fitted in order to match the observations of the positron 




FIGURE 2. Positron fraction from the decay of gravitinos 
for a NFW halo profile and the M2, MED and Ml propagation 
models. Gravitino parameters are as in Fig.[T] We also show for 
comparison the contribution to the positron fraction from spal- 
lation of cosmic rays on the Galactic disk, which constitutes the 
background to our signal. 



fraction. For the flux of secondary positrons, and the pri- 
mary and secondary electrons we used the parametriza- 
tions obtained in |31] from detailed computer simula- 
tions of cosmic ray propagation |33]. 

We show in Fig. |2]the positron fraction for different 
diffusion models in the case of the NFW profile when 
m3/2 — 150GeV and T3/2 — 1.3 x lO^^s. Interestingly, 
we find that gravitino parameters which predict a de- 
parture from a simple power law in the extragalactic 
gamma-ray spectrum at energies above 2 GeV (as ob- 
served by EGRET), inevitably predict a bump in the 
positron fraction at energies above 7 GeV (as observed 
by HEAT) |17fl Furthermore, the presence of this fea- 
ture is not very sensitive to the many astrophysical un- 
certainties. This remarkable result holds not only for the 
scenario of gravitino dark matter with broken /^-parity, 
but also for any other scenario of decaying dark matter 
with lifetime ~ 10^^ s which decays predominantly into 
Z" and/or gauge bosons with momentum ~ 50 GeV. 

Antiproton flux 

The general antimatter transport equation can be sim- 
plified by taking into account that energy losses are negli- 
gible for antiprotons. Using this approximation, we have 
solved the diffusion equation for three propagation mod- 
els consistent with the B/C ratio which yield the mini- 
mum (MIN), median (MED) and maximal (MAX) an- 
tiproton flux 1I30I1 . We have found that the Green's func- 
tion which describes antiproton propagation from dark 
matter decay can be approximated by the following in- 



' The same conclusion has been independently reached by Ishiwata, 
Matsumoto and Moroi in fT^ . 



TABLE 2. Coefficients of the interpolating 
function Eq. l[8j for the antiproton Green's 
function assuming a NFW halo profile and for 
the MIN, MED and MAX propagation models 
defined in | 30]. 



model 


X 


y 


z 


MIN 


-0.0537 


0.7052 


-0.1840 


MED 


1.8002 


0.4099 


-0.1343 


MAX 


3.3602 


-0.1438 


-0.0403 



terpolating function 111711 : 

Gp{Tj') - 10^4e-^+-^''"^+-''"'^5(r'-r)cm-3s , (8) 

which, again, is valid for any decaying dark matter par- 
ticle. The coefficients x, y and z for the NFW profile can 
be found in Table|2]for the different diffusion models. In 
this case the approximation is accurate to a 5-10%. As in 
the case of the positrons, the dependence of the Green's 
function on the halo model is fairly weak. 

The interstellar antiproton flux can be then straight- 
forwardly computed from Eqs. (|4|l and Eq. (|5]l using 
the previous parametrization of the Green's function and 
the energy spectrum of antiprotons from gravitino decay, 
dNp/dT. However, this is not the antiproton flux mea- 
sured by balloon or satellite experiments, which is af- 
fected by solar modulation. In the force field approxima- 
tion 13411 the effect of solar modulation can be included 
by applying the following simple formula that relates the 
antiproton flux at the top of the Earth's atmosphere and 
the interstellar antiproton flux ll35ll : 



<i>r(7TOA) 



ImpT-YOA + T^ 



TOA 



2m„Tis 



^IS 



^'p{Tis), (9) 



where Tjs — Tjo/^ + (j)p, with Tjs and Txoa being the 
antiproton kinetic energies at the heliospheric boundary 
and at the top of the Earth's atmosphere, respectively, 
and (j)F the solar modulation parameter, which we take 
(pF = 500 MV. 

We show in Fig. [3] the predicted antiproton flux from 
gravitino decay for the MIN, MED and MAX diffusion 
models when 1112^2 — 150GeV and T3/2 — 1-3 x lO^^s. 
From the plot, the extreme sensitivity of the primary an- 
tiproton flux to the choice of the diffusion model is appar- 
ent: parameters that successfully reproduce the observed 
B/C ratio lead to antiproton fluxes that span over two or- 
ders of magnitude. For a wide range of propagation pa- 
rameters, the total antiproton flux is well above the ob- 
servations and thus our scenario is most likely excluded. 
However, the MIN model yields a primary flux that is be- 
low the measured flux and thus might be compatible with 
observations. 

We have analyzed more carefully the predictions for 
the MIN model computing the different contributions to 
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FIGURE 3. Contribution to the antiproton flux at the top of 
the atmosphere from the decay of dark matter gravitinos for the 
NFW halo profile and the MIN, MED and MAX propagation 
models. Gravitino parameters are as in Fig.[T] 



the total antiproton flux. The result is shown in the right 
panel of Fig |4] where, for consistency, we have adopted 
as background the secondary antiproton flux calculated 
in 113611 for the same MIN model. Although the primary 
antiproton flux is smaller than the measured one, the to- 
tal antiproton flux is a factor of two above the observa- 
tions. Nevertheless, in view of all the uncertainties that 
enter in the calculation of the antiproton flux, it might 
be premature to conclusively rule out the scenario of de- 
caying gravitino dark matter. Namely, in addition to the 
uncertainties stemming from degeneracies in the diffu- 
sion parameters, there are also uncertainties from the nu- 
clear cross sections and, to a lesser extent, uncertainties 
from the description of the interstellar medium and solar 
modulation (for a discussion of the various uncertainties 
see Ii36]). Furthermore, we used a simplified diffusion 
model that neglects the effects of reacceleration, energy 
losses and tertiary contributions. Therefore, there could 
be certain choices of parameters or more refined diffu- 
sion models where the total antiproton flux is consistent 
with experiment^. 



CONCLUSIONS 

The scenario of gravitino dark matter with broken R- 
parity naturally reconciles three popular paradigms that 
seem to be in mutual conflict: supersymmetric dark mat- 
ter. Big Bang Nucleosynthesis and thermal leptogenesis. 
Moreover, the gravitino decay products might become 
observable, thus opening the possibility of the indirect 
detection of the elusive gravitino dark matter. 



^ Some works have reported a deficit in the predicted secondaiy an- 
tiprotons compai'ed to the observations and argued that this deficit 
could be connected with a contribution of primary antiprotons f37ll . 
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FIGURE 4. Summary of the signatures of gravitino dark matter decay in the extragalactic gamma-ray flux (left), the positron 
fraction (middle) and the antiproton flux (right), compared to the EGRET, HEAT and BESS data respectively. The contribution 
from gravitino decay is shown with a dotted line, the brackground with a dashed line and the total flux with a solid line. In these 
plots, we have adopted the MIN diffusion model, m-^j2 ^ ISOGeV and T3/2 ~ 1.3 x 10^^ s. 



We have shown that the EGRET anomaly in the ex- 
tragalactic gamma-ray flux and the HEAT excess in the 
positron fraction can be simultaneously explained by the 
decay of a gravitino with a mass m3/2 ^ 150 GeV and 
a lifetime of T3/2 ~ 10^^ s. However, the predicted an- 
tiproton flux tends to be too large, although the prediction 
suffers from large uncertainties and might be compatible 
with present observations for certain choices of propaga- 
tion parameters. Our results are summarized in Fig.|4] 

This remarkable result holds not only for the scenario 
of gravitino dark matter with broken 7? -parity, but also for 
any other scenario of decaying dark matter with lifetime 
~ 10^^ s which decays predominantly into Z" and/or 
gauge bosons with momentum ~ 50 GeV. 
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